arXivil 508.0547 lvl [astro-ph.GA] 22 Aug 2015 


Publ. Astron. Soc. Japan (2014) 00(0), 1-16 
doi: 10.1093/pasj/xxx000 


1 


Constraints on the Minor Merging and Star 
Formation History of the Wolf-Rayet Galaxy 
NGC 5430 Through Observations 

Hsi-An Pan, Nario Kuno, Kazuo Sorai and Michiko Umei 

Department of Physics, Faculty of Science, Hokkaido University, Kita 10 Nishi 8 Kita-ku, 
Sapporo 060-0810, Japan 

Faculty of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennoudai, Tsukuba, 
Ibaraki 350-8577, Japan 

Department of Physics, Faculty of Science, Hokkaido University, Kita 10 Nishi 8 Kita-ku, 
Sapporo 060-0810, Japan 

Department of Physics, Faculty of Science, Hokkaido University, Kita 10 Nishi 8 Kita-ku, 
Sapporo 060-0810, Japan 

*E-mail: hapan@astro1 .sci.hokudai.ac.jp 
Received ; Accepted 

Abstract 

We used multi-wavelength analysis of the newly observed molecular gas ( 12 CO and 13 CO (1- 
0)) with interferometer CARMA and archival star formation tracers to constrain the interaction, 
merging, and star formation history of an off-center minor merger, a three-spiral barred galaxy 
NGC 5430 and its satellite embedded in the bar. Morphology of the molecular gas in the 
bar of NGC 5430 shows minimal signs of recent interactions in our resolution. The apparent 
morphological remnant of the past galaxy interaction is an asymmetric spiral arm, containing 
more molecular gas and exhibiting higher star formation rate (SFR) surface density than the 
two primary arms. Rotation curve analysis suggests that NGC 5430 and its satellite collided 
several Gyr ago. History of star formation was constrained by using SFRs that trace differ¬ 
ent timescales (infrared, radio continuum, and Ha). The collision occurred 5-10 Myr ago, 
triggering a transient off-center starburst of Wolf-Rayet stars at the eastern bar end. In the 
past, the global SFR during the Wolf-Rayet starburst peaked at 35 M 0 yr _1 . At present, the 
merger-driven starburst is rapidly decaying and the current global SFR has decreased to the 
Galactic value. The SFR will continue to decay as suggested by the present amount of dense 
gas (traced by HCN (1-0)). Nonetheless, the global SFR is still dominated by the Wolf-Rayet 
region rather than the circumnuclear region. Compared with other barred galaxies, the circum- 
nuclear region exhibits particularly low dense gas fraction, low star formation activity and high 
concentration of gas. Physical properties of the molecular gas are inferred by using the large 
velocity gradient (LVG) calculations. The initial mass ratio of the NGC 5430 and its satellite are 
suggested to be in the intermediate ratio range of 7:1-20:1. 

Key words: Galaxies: individual: NGC 5430 - Galaxies: interactions - ISM: molecules - Galaxies: star 
formation 
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1 Introduction 

Hierarchical formation models suggest that galaxy-galaxy in¬ 
teractions play a key role in the formation, evolution, and mor¬ 
phology of galaxies (e.g. White & Rees 1978; Kazantzidis et 
al. 2008). Minor mergers, i.e., mergers of galaxies with large 
mass ratio, take place frequently in z < 1, and their rate is sig¬ 
nificantly higher than that of major mergers (mergers between 
galaxies of equal mass). While major mergers tend to form fea¬ 
tureless massive elliptical galaxies, minor mergers are more rel¬ 
evant to disk galaxies, because they exert a relatively small ef¬ 
fect on the primary disks (e.g. Mastropietro et al. 2005; Casteels 
et al. 2014). Minor mergers help to form bulge and increase 
disk thickness in disk galaxies through deposition of satellite 
materials (Shioya et al. 1998; Aguerri et al. 2001; Eliche-Moral 
et al. 2006; Hopkins et al. 2010). These mergers also help to 
form galactic bars by triggering disk instability (Hernquist & 
Mihos 1995; Skibba et al. 2012), and probably trigger some 
low or intermediate luminosity active galactic nuclei (AGNs) 
(Hopkins et al. 2008; Kazantzidis et al. 2008; Villalobos & 
Helmi 2008; Mapelli et al. 2015). 

Galaxy-galaxy interactions and mergers trigger star forma¬ 
tion. Observations suggest that galaxy-galaxy interactions in¬ 
creases star formation activity through violent compression of 
molecular gas (Lambas et al. 2003; Kauffmann et al. 2004; Lin 
et al. 2007; Smith et al. 2007; Woods & Geller 2007; Kaviraj 
2014). Kaviraj (2014) suggested that about half of the star for¬ 
mation events in the local Universe are triggered by the minor 
merger processes. The notion of merger-induced star forma¬ 
tion is also supported by theoretical studies (e.g. Cox et al. 
2008; Lotz et al. 2010). 

Observations suggest that off-center mergers are common 
and can occur repeatedly in a primary disk (e.g. Mandelker 
et al. 2014; Guo et al. 2015). The collision sites are shown as 
discrete star forming regions in the disks of primary galaxies 
and contribute significant amount of UV and optical light to the 
host galaxies (Wuyts et al. 2012; Guo et al. 2015). 

In this work, we diagnosed the history of mergers and star 
formation events in a nearby off-center minor-merger system: 
NGC 5430 and its satellite galaxy HOLM 569B 1 . The results 
can help to constrain the star formation properties and history 
of the high-redshift counterparts that often exhibit higher rate of 
mergers but cannot be resolved in details. 

A Sloan digital sky survey (SDSS) g — r — i composite image 
of the primary disk galaxy NGC 5430 is shown in Ligure 1. 
NGC 5430 (14h00m45.7s, +59° 19'42") is classified as SB(s)b, 
located at 42 Mpc, where l" = 205 pc. NGC 5430 is a member 
of a galaxy group (Geller & Huchra 1983). The satellite dwarf 
galaxy HOLM 569B (14h00m47.3s, +59° 19 / 27 // ) is embedded 
in the eastern bar end, and is shown as a blue spot in Ligure 1, 
around 4.5 kpc (22") from the galaxy center (Keel 1982). A 2D 

1 In the Holmberg galaxy catalog of galaxy groups, NGC 5430 is HOLM569. 


Gaussian fit to its i-band suggests that the diameter of HOLM 
569B is ~ 1 kpc. 

The minor-merger system is manifested by the asymmetric 
spiral arms and the unusual star formation activity at the eastern 
bar end. NGC 5340 contains three spiral arms, including two 
primary arms and one asymmetric arm located at the north of 
the galactic nucleus, implying a morphology distortion caused 
by an external effect. The global infrared luminosity of NGC 
5430 is not particularly high (~ 7.5 x 10 10 L©; Sanders et al. 
2003) and is inferior to the luminous infrared galaxy (LIRG; > 
10 11 L©). However, the giant HII region at the eastern bar end 
is unusually bright, with star formation rate equal to that of the 
entire Milky Way. Moreover, the giant HII region deviates from 
the global HII regions luminosity function of the galaxy, sug¬ 
gesting that the bright core is a separate object (Keel 1982). The 
most massive stars are forming in the giant HII region, includ¬ 
ing 10 4 -10 5 OB stars and Wolf-Rayet stars (> 20 M©) (Keel 
1982; Lernandes et al. 2004). In the rest of this work, we re¬ 
fer to the eastern bar end (the position of the HOLM 569B and 
the surrounding 1-kpc-diameter area) as the Wolf-Rayet (W-R) 
region. The contribution of UV light and the size of the W-R 
region are similar to those of the giant star forming clumps that 
are formed through minor mergers in the local and high-redshift 
Universe (Elmegreen et al. 2007; Livermore et al. 2012; Guo et 
al. 2015). 

The spatial distribution of molecular gas in the NGC 5430 
was unknown prior to this work, while single-point observations 
were performed in a few studies. Kruegel et al. (1990) observed 
low density gas in 12 CO (1-0), 12 CO (2-1) and 13 CO (2-1) in 
the galactic nucleus by using the IRAM 30-m telescope. Later, 
Contini et al. (1997) used the same telescope to measure dense 
gas tracers HCN (1-0) and CS (3-2), as well as to perform 12 CO 
(2-1) observations towards the nucleus and the W-R region. No 
CS (3-2) emission was detected in both regions, while HCN (1- 
0) was detected in the nucleus. The HCN-far infrared (LIR) lu¬ 
minosity relation for this nucleus is the same as for other nearby 
galaxies; however, the 12 CO (2-1), as well as the 12 CO (1-0) 
intensities are unusually high relative to the LIR fluxes. The 
unusual nuclear CO intensities were reported by Contini et al. 
(1997) for two other Wolf-Rayet galaxies. 12 CO (2-1) was de¬ 
tected in the W-R region. The intensity of 12 CO (2-1) in the 
W-R region was 13 times weaker than in the nucleus. 

This paper is structured as follows. The mapping observa¬ 
tions of 12 CO (1-0) and 13 CO (1-0) are described in Section 
2. The results of these observations, including the images, mor¬ 
phological and kinematic analyses, are shown in Section 3. In 
Section 4 , we present the physical properties of molecular gas 
and the history of star formation. In Section 5, we discuss the 
possible initial mass ratio of NGC 5430 and its satellite galaxy. 
In Section 6, we summarize our work. 
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Fig. 1 . SDSS (Baillard et al. 2011) g — r — i composite of NGC 5430 system 
taken from the NED. The white circle indicates the area of CARMA observa¬ 
tions. Diameter of the circle is 80" or ~ 16 kpc. 


2 Observation and Data Reduction 

Molecular gas of NGC 5430 was mapped by using the 
combined array for research in millimeter-wave astronomy 
(CARMA). CARMA consists of six 10.4-m-diameter anten¬ 
nas and nine 6.1-m-diameter antennas. Observations were per¬ 
formed on July 30th, August 1st and 3rd of 2013, by us¬ 
ing E configuration, with baselines in the 8 - 66 m range. 
Observations were performed towards the central 80" (diam¬ 
eter of ~ 16 kpc) of the galaxy (white circle in Figure 1). Two 
narrow bands (250 MHz) were centered at 115.27 GHz ( 12 CO 
(1-0)) and 110.20 GHz ( 13 CO (1-0)). The total velocity width 
of a narrow band was ~ 740 km -1 , with spectral resolution of 
~ 9 km s -1 . Wide bands (500 MHz) were placed between the 
narrow bands for the calibration of continuum. Flux, passband, 
and phase calibrators were mwc349, 1635+381, and 1642+689, 
respectively. The total observation time was ~ 15 hours, includ¬ 
ing the target and the calibrators. 

Calibration, imaging and deconvolution were performed by 
using the standard procedures of the MIRIAD package (Sault 
et al. 1995). The final 12 CO (1-0) and 13 CO (1-0) clean maps 
have the spatial resolutions of 7.6" (1.5 kpc) x 5.3" (1.1 kpc) 
with RA. = -71.8°, and 8.1" (1.7 kpc) x 5.5" (1.1 kpc) with 
P.A. = -68.0°, respectively. With the velocity resolution binned 
to 12 km s -1 , the final spectral cube of 12 CO (1-0) has the 
sensitivity of ctrms ~ 11 mJy beam -1 . Because 13 CO (1-0) 
is significantly weaker than 12 CO (1-0), it has been binned to 
the lower spectral resolution of 25 km s -1 for improving the 
signal-to-noise ratio (S/N). The final ctrms of 13 CO (1-0) is 6 
mJy beam -1 . The largest structure that can be detected in these 
observations is ~ 82" or ~ 17 kpc. Therefore missing flux is 
negligible for this observation. 


3 Results 

3.1 Galactic Morphology 

3.1.1 Channel Maps of CO 

The channel map of 12 CO (1-0) (hereafter 12 CO) is delineated 
in Figure 2 with white contours. The SDSS i-band image is 
also shown, in both greyscale and black contours. The galac¬ 
tic center and the W-R region are marked with crosses, respec¬ 
tively. The 12 CO exhibits a wide range of velocities, from 2730 
to 3170 km s -1 . The strongest emission was generated at the 
galactic center (~ 40 cr RM s or ~ 0.44 Jy beam -1 ). Two elon¬ 
gated structures emerge from the center, and were detected with 
the significance of ~ 10 ctrms or ~ 0.11 Jy beam -1 . The elon¬ 
gated structures correspond to the bar regions of the optical im¬ 
age. The W-R region was detected with the significance of ~ 
10 ctrms or ~ 0.11 Jy beam -1 . 

The channel map of 12 CO reveals features that are typical 
of a barred galaxy. The 12 CO emission from the bar shifts to¬ 
wards the downstream of galactic rotation, known as the offset 
ridges of the bar. The offset ridges trace the shocks that are 
associated with the crowded gas orbits, as suggested by theoret¬ 
ical studies (Athanassoula 1992). The two offset ridges appear 
to be spatially and kinematically symmetric with respect to the 
galactic center. After passing through the offset ridge shocks, 
the gas moves inwards and accumulates in the central region, 
yielding high gas concentration at the center (e.g., Kuno et al. 
2007). This phenomenon is well-known as bar-driven gas in¬ 
flow (Sakamoto et al. 1999). The widths of velocity ranges of 
the galactic center, eastern, and western offset ridges were ~ 
440, 180, and 160 km s -1 , respectively. 12 CO was also de¬ 
tected in the asymmetric arm, emerging in the 2780 to 2860 
km s -1 range. 12 CO was observed in a spot of the north pri¬ 
mary arm with velocity in the 2854 - 2830 km s -1 range and 
with significance of ~ 4 ctrms or 0.044 Jy beam -1 . 

Figure 3 shows the channel map of 13 CO (1-0) (hereafter 
13 CO). The keys are the same as in Figure 2. 13 CO emission 
was generated in the same velocity range as that of 12 CO. The 
galactic center was solidly detected (~ 14 cr RM s or 0.084 Jy 
beam -1 ). The eastern bar was detected with the significance of 
~ 8 ctrms or 0.048 Jy beam -1 , while the other side was barely 
detected with ~ 3 ctrms or 0.018 Jy beam -1 . The 13 CO of 
the asymmetric arm was observed in the velocity of ~ 2800 km 
s -1 , consistent with the velocity of 12 CO emission. We did not 
detect any 13 CO from the two primary arms. 

3.1.2 Integrated Intensity Maps of CO 

Integrated intensity maps of 12 CO and 13 CO are shown in 
Figure 4(a) and 4(b) with white contours, overlaid on the SDSS 
i-band image. The morphology of molecular gas at the center 
and throughout the bar appears to be not severely disturbed by 
the galactic mergers. Both 12 CO and 13 CO intensity maps con¬ 
tain a bright central component with the radius of ~ 1.5 kpc. 
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Fig. 2. Channel map of 12 CO (1-0) in white contours overlaid on the SDSS i-band image (Baillard et al. 2011) in black contours and grey scale. Velocity 
resolution of the data cube is 12 km s -1 . The figure presents every other channel to save the space. The contours of 12 CO are 2, 3, 4, 5, 7, 10, 15, 20, and 30 
cr rms, where 1 ctrms is 11 mJy. The crosses denote the galactic center (upper one) and the W-R region (lower one), respectively. Beam size of 12 CO (1-0) 
is 7 . 6 " x 5.3" (1.5 x 1.1 kpc) with RA. = -71.8°, showing in the lower-right corner of each channel. Velocity of each channel is displayed at the upper-left. 


There is an offset between the strongest CO intensity and the 
galactic center. The reason for this is enclosure of the galactic 
nucleus within a ring-like structure (Garcia-Barreto et al. 1996); 
thus, the CO peak intensity position is determined from the con¬ 
volution of the ring with our CO beam. Therefore, the CO peak 
is not necessarily at the galactic center if the gas distribution in 
the ring is not uniform. Offset ridges at the two sides of the 
bar are seen in 12 CO, and are symmetric with respect to the 
galactic center. The offset ridges extend from the center to the 
distance of ~ 5.6 kpc, and have the width of ~ 1.2 kpc. The 
eastern ridge exhibits stronger 12 CO intensity than the western 
ridge. Within this resolution, it cannot be determined whether 
the eastern ridge is, in general, stronger throughout the entire 
ridge or whether the W-R region, containing extra molecular 
gas contributed by the satellite, has been convolved with a large 
beam. 13 CO is only detected in the eastern ridge. The length 
and width are consistent with those observed for 12 CO. Because 


the two sides of the bar exhibit different star formation activity 
(stronger at the eastern ridge and weaker at the western ridge), 
their intrinsic 12 C, 13 C, 12 CO, and 13 CO abundances can be 
different. Both 12 C and 13 C can be enriched by the rapid cy¬ 
cling of interstellar medium (ISM) through stars. Once they are 
used to form 12 CO and 13 CO, the molecules can be destroyed 
by photo-dissociation through the UV photons from stars. The 
intrinsic 12 CO/ 13 CO abundance ratio then depends on the net 
result of these processes (e.g., Milam et al. 2005; Henkel et al. 
2014; Szucs et al. 2014 and references therein). The observed 
line intensity ratio of 12 CO/ 13 CO also depends on the radiative 
transfer and on the excitation of molecular lines. Therefore, 
star formation activity and the ISM environment play critical 
roles in determining the observed 12 CO/ 13 CO ratio. The unde¬ 
tected 13 CO of the western ridge may imply non-similar physi¬ 
cal properties of molecular gas across the two sides of the bar. 


The asymmetric spiral arm is connected to the eastern ridge 
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Fig. 3. Channel map of 13 CO (1-0) in white contours overlaid on the SDSS i-band image in black contours and grey scale. The symbols are the same as in 
Figure 2. The contours are 2, 3, 4, 5, 7, and 10 ctrms, where 1 o-rms is 6 mJy. Beam size of 13 CO (1-0) is 8.1" x 5.5" (1.7 xl.1 kpc) with P.A. = -68.0°. 
Note that the channel velocity width of 13 CO (1-0) is larger than 12 CO (1-0). Velocity resolution of 13 CO (1-0) data cube is 25 km s _1 . 


in 12 CO. 12 CO emission of the asymmetric spiral arm is spa¬ 
tially correlated with the stellar light traced in SDSS image. The 
asymmetric spiral arm is more fragmented in 13 CO, and is not 
connected to the eastern ridge. The detected 13 CO is found in 
the regions with high 12 CO intensity and stellar light. 

3.1.3 Central Concentration of Molecular Gas and Bar 
Strength 

The total gas mass (H 2 + He) in the central ~ 1 kpc of NGC 
5430 is considerably higher than other strong bar (SB and SAB) 
galaxies. We estimated the mass of gas in the central ~ 1 kpc (in 
diameter) using the flux of the beam (~ 110 Jy beam -1 km s -1 ) 
at the galactic center, along with the Galactic CO-to-H 2 conver¬ 
sion factor (Xco) of 2 x 10 20 cm -2 (K km s -1 ) -1 (Bolatto 
et al. 2013), and a correction factor of 1.36 for He and heavy 
elements. Here we assumed that Xco is constant and similar to 
that of the Milky Way disk across entire NGC 5430 when the 
gas properties are averaged over several 100 parsec to kilopar- 
sec scale (e.g., Donovan Meyer et al. 2013; Bolatto et al. 2013). 
The derived total mass is ~ (8.3 ± 1.5) x 10 s M© 2 . The value 
is considerably higher than the average value of galaxies with 
the same morphological type, ~ 3.3 x 10 8 M 0 (Sakamoto et 

2 The uncertainty is dominated by the uncertainty of flux calibration of 
CARMA observations. 


al. 1999) 3 . 

Kuno et al. (2007) show that the central concentration of 
molecular gas is proportional to the bar strength because strong 
bar exerts large gravitational torque and be most effective in 
driving gas inflows (see their Figure 49). Central concentration 
of molecular gas (/ co n) is defined as: 

j. _ Mh 2 (72k2o/8) ^ 

/con_ M H2 (fW2)’ W 

where Rk 2 o is the radius of 20 mag/arcsec -2 in the 2MASS 
Tfs-band, and Mh 2 (-RK 20 / 8 ) and Mh 2 ( 77 k 2 o/ 2 ) represent the 
gas mass enclosed within the radius of 77 k 2 o /8 and f?K 2 o/ 2 , re¬ 
spectively. Using Equation (1), the f con of NGC 5430 is esti¬ 
mated as ~ 0.34. The value of f con is indeed higher that the 
average value (0.20) of barred galaxies in Kuno et al. (2007) 4 . 
Therefore, NGC 5430 likely has a stronger bar than the visually 
isolated galaxies in Kuno et al. (2007) in terms of f con . 

The strong bar is also implied by the high de-projected bar 
ellipticity. There is increasing evidence that bar morphology 
is determined primarily by the bar strength. Various studies, 
including both observations and simulations, have reported the 

3 Corrected for the adopted Xco in this work. 

4 Galaxies in Kuno et al. (2007) have distance of < 25 Mpc and an observing 
resolution of 15". NGC 5430 is farther than those galaxies by ~ 2 times, 
but observed with a resolution ~ 2 times finer. Thus these two studies 
sample the same galactocentric ranges in Equation (1). 
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(a) 



14 h 0 m 50 s 48 s 46 s 44 s 42 s 


(b) 

Fig. 4. (a) Integrated intensity map of 12 CO (1-0) (white contours) overlaid 
on SDSS i-band image (black contours and grey scale). Contour of 12 CO 
(1-0) is plotted with steps of 0.6, 0.9, 1.5, 1.9, 3, 4, 6, 8, 13, 20, 30, 40, 
60, 80, 100 Jy beam -1 km s -1 . The galactic center and the W-R region 
are marked with crosses. The noise level of the intensity map is ~ 0.2 Jy 
beam -1 km s -1 . (b) Integrated intensity map of 13 CO (1-0) (white con¬ 
tours) superimposed on SDSS i-band image. Contour of 13 CO is displayed 
with steps of 0.3, 0.4, 0.9, 1.7, 2.0, 3.0, 4.5, 6.0, 8.0, 10, 13 Jy beam -1 km 
s -1 . The noise level is ~ 0.1 Jy beam -1 km s -1 . 

existence of correlation between the bar strength and bar ellip- 
ticity (Athanassoula 1992; Martinet & Friedli 1997; Abraham & 
Merrifield 2000; Block et al. 2001; Block et al. 2004; Comeron 
et al. 2010; Kim et al. 2012). The Spitzer project of S 4 G (Sheth 
et al. 2010) has decomposed the galactic structure of NGC 5430 
by using the GALFIT package and the Spitzer IRAC 3.6 pm 
image, by assuming a three-component bulge-disk-bar model 
(Salo et al. 2015). The derived bar ellipticity is 0.6. The bar 
ellipticity of NGC 5430 is again on the higher side of the ob¬ 
served and simulated bar ellipticity, which ranges from 0.2 to 
0 . 8 . 

It is certainly possible that the formation of the strong 


bar in NGC 5430 is induced by the past galaxy interaction. 
Simulations have shown that even a small perturbation, such as 
galaxy flyby and minor merger, can induce a bar (e.g., Noguchi 
1987; Berentzen et al. 2004; Lang et al. 2014; Lokas et al. 
2014). 

3.2 Velocity Fields 

3.2.1 Velocity Fields of CO and the Comparison with Ha 

Figures 5(a) and 5(b) show the velocity fields of 12 CO and 
13 CO, respectively. Two crosses mark the galactic center (upper 
cross) and the W-R region (lower cross). The receding and ap¬ 
proaching sides are located in the south and north, respectively. 
The velocity fields of 12 CO and 13 CO are consistent with each 
other on the global scale. 

The detailed velocity distribution in the 13 CO field is af¬ 
fected by the low sensitivity. For this reason, the kinematic dis¬ 
cussion will focus on the 12 CO field. The central region reveals 
parallel velocity contours in 12 CO. However, velocity fields of 
disk galaxies usually exhibit either the “spider” diagram of cir¬ 
cular motion or the “S-shape” twisted velocity gradient of non¬ 
circular motion. Our CO emission may be smeared out owing to 
the large beam size and low velocity resolution (12 km s _1 ) at 
the galactic center, where velocity field changes rapidly within 
a single beam. Therefore neither the spider nor the S-shape pat¬ 
tern is seen. It is marginally observed from the 12 CO map that 
the velocity gradient is almost normal to the CO offset ridges 
(namely the contours are parallel to the ridges), implying a large 
velocity jump across the bar, which is predicted by theoretical 
studies (e.g., Athanassoula 1992). For 12 CO, this velocity jump 
is about 100 km s -1 per 1 kpc. 

Figure 5(c) compares the velocity field of 12 CO (white con¬ 
tours) to the Ha from the GHASP survey (color scale and black 
contours). Spatial and velocity resolutions of the Ha field are 
~ 4" (820 pc) and 5 km s _1 , respectively. The CO field reveals 
the global variation in velocity that is observed in Ha. Because 
of the higher resolution, the velocity field of Ha reveals more 
velocity structures at the galactic center. The galactic center re¬ 
veals twisted S-shape velocity contours, which is an evidence of 
the bar non-circular motion (Epinat et al. 2008). This pattern is 
barely seen if the velocity field of Ha is smoothed to the veloc¬ 
ity and spatial resolutions of the CO maps (Figure 5(d)). Figure 
5(e) shows the velocity residual (color scale) of the smoothed 
Ha field with the CO field subtracted from it, overlaid on the 
emission boundary of the 12 CO integrated intensity map. The 
majority of the residuals are within ±10 km s -1 , while some 
regions deviate by as much as ± 50 km s -1 . The regions with 
large residuals are located around the edge of the CO emissions 
or are associated with weak CO emission. Thus, these large 
residuals are not significant for judging the difference between 
molecular and ionized gases. 
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Velocity dispersions of molecular and ionized gases appear 
to be more different. Figure 6(a) shows the velocity dispersion 
of the molecular gas, while panels (b) and (c) show the velocity 
dispersion of the ionized gas in color scales. Velocity dispersion 
of the molecular gas reaches the peak of ~ 90 km s _1 at the 
galactic center, with an average of ~ 65 km s -1 within the cen¬ 
tral 1 kpc. In the bar region, large velocity dispersion appears 
at the leading side of the bar. The eastern bar exhibits a slightly 
larger dispersion than the western bar, ~ 15-40 km s -1 com¬ 
pared with ~ 12-30 km s _1 . The asymmetric spiral shows a 
low velocity dispersion of 5 - 12 km s -1 . We found that this dis¬ 
tribution of velocity dispersion is consistent with that observed 
in low-resolution 12 CO and can be explained by secular origin. 
The high velocity dispersion (> 50 km s -1 ) at the galactic cen¬ 
ter results from the accumulation of molecular gas in the colli¬ 
sion region where x\ (bar) and X 2 (circumnuclear region) orbits 
intersect, and from the spatial superposition of circular and non¬ 
circular components (e.g., Huttemeister et al. 2000; Schinnerer 
et al. 2000; Hsieh et al. 2011). The intermediate velocity dis¬ 
persion (a few tenths km s -1 ) of the bar, peaking at the leading 
side on both sides of the bar, is owing to the crowded obits (e.g., 
Athanassoula 1992; Huttemeister et al. 2000). The low velocity 
dispersion (< 15 km s -1 ) in the asymmetric spiral arm is con¬ 
sistent with that of the unresolved giant molecular clouds (asso¬ 
ciations) in the galactic disks (e.g., Tosaki et al. 2007; Donovan 
Meyer et al. 2013; Faesi et al. 2014). 


In contrast to the CO, velocity dispersion of the ionized gas 
exhibits stronger local variation because of the higher resolu¬ 
tion. Figures 6(b) and 6(c) show the velocity dispersion of Ha 
using color scale. The CO boundary is overlaid with red con¬ 
tour in Figures 6(b). Two regions with high velocity disper¬ 
sion are observed in the central region, located at the opposite 
sides of the galactic center. These regions correspond to the 
intersection of x\ and X 2 orbits, the so-called contact points 
(indicated by two white lines in Figure6(c)). The velocity dis¬ 
persion of Ha peaks at the northern intersection with 72 km 
s- 1 , while the maximal velocity dispersion at the southern in¬ 
tersection is 64 km s -1 . The intersections are connected by two 
high-dispersion dust lanes of the bar. The maximal velocity dis¬ 
persion at the dust lanes reaches ~ 65 km s -1 . Beyond the dust 
lanes, the W-R region exhibits larger velocity dispersion than 
the opposite region in the western bar. The peak velocity dis¬ 
persion in the W-R region is 61 km s -1 with an average of 45 
km s -1 within the surrounding 1 kpc region. The large velocity 
dispersion of Ha compared with the CO is perhaps owing to 
the bias towards the diffuse ionized gas disturbed by supernova 
explosions and the expanding nebulae around the massive pop¬ 
ulations, such as the W-R stars (Melnick et al. 1999; Moiseev & 
Lozinskaya 2012; Kam et al. 2015). 


3.2.2 Position-Velocity Diagrams and Rotation Curve 

Figure 7 compares the position-velocity (PV) diagrams of Ha 
(color scales, grey contours), and 12 CO (black contours). The 
PV diagrams are obtained along the major axis of the velocity 
fields at 181.9°. The PV diagram of the ionized gas has the 
same resolution as that of the molecular gas. Both PV diagrams 
are composed of a circumnuclear disk with steep increase in 
velocity within d= 20" (~ 4 kpc), and a strong emission at ~ 
30" 6 kpc) and 2800 km s -1 tracing the asymmetric spiral 

arm. 

The blue curve in Figure 7 represents the rotation curve of 
the NGC 5430 measured in the ionized gas by Epinat et al. 
(2008). For comparison between the receding and the approach¬ 
ing sides, the reversed rotation curve is overlaid with pink dots. 
The rotation curve reaches a plateau within a few arcsec on both 
sides, fluctuating at about 3100 km s -1 (+ 140 km s -1 relative 
to the V sys of 2960 km s -1 ) and 2800 km s -1 (- 160 km s -1 ). 
The two sides of the rotation curve are symmetric out to ~ d= 

U 

20 (rsj 4 kpc). This radius includes the central region and ma¬ 
jority of the bar. A small bump is observed at the receding side 
at the radius just beyond the W-R region, ~ 30" (~ 6 kpc), 
where the velocity at the receding side is ~ 10 to 20 km s -1 
larger than that at the approaching side. The kinematic asym¬ 
metry around the W-R region might be owing to the inhomo¬ 
geneous mass distribution caused by extra gas contributed by 
the satellite collisions and the subsequent active star formation 
(note that the rotation curve is derived from Ha and therefore it 
is sensitive to the star formation history). 

In terms of the rotation curve, the NGC 5430 has almost 
relaxed and shows nearly no signs of interaction in the inner 
disk. Overall, we found that the rotation curve beyond the cen¬ 
tral region (> |=blkpc|) is flat, fluctuating within ~ 60 kms -1 
on the both sides. Moreover, the extent of asymmetry between 
the two sides is small (< 50 km s -1 ). Simulations reveal that 
systems with ongoing minor galaxy interaction or encounters 
have non-flat rotation curves. Velocity increases with galac- 
tocentric radius by at least 100 km s -1 in the galactic disk. 
Moreover, the rotation curves are asymmetric, differing by > 
100 km s -1 between the two sides of the disk (Kronberger et al. 
2006). Therefore, the observed rotation curve of NGC 5430 im¬ 
plies that the system has almost relaxed, exhibiting no obvious 
signs of kinematic distortion. 

The outer disk shows a larger degree of asymmetry than the 
inner disk. Beyond the asymmetric spiral arm (d= 40" or ~ 8 
kpc), the two sides of the rotation curve exhibit a larger devia¬ 
tion by ~ 30 km s -1 . It is possible that the outer disk is not com¬ 
pletely stabilized from the past interaction with HOLM569B. 
Moreover, because NGC 5430 is the member of a galaxy group 
(Geller & Huchra 1983), the perturbed outer disk can also result 
from ongoing interaction with other group members. 

The kinematics of NGC 5430 helps to constrain the 
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(a) (b) (c) 

Fig. 6. Velocity dispersions of NGC 5430. The galactic center and the W-R 
region are marked with black crosses in all panels (a) Velocity dispersion 
(second moment) of 12 CO (1-0). (b) Velocity dispersion of Ha (color scale) 
with the boundary of 12 CO integrated intensity map overlaid (red contour), 
(c) Velocity dispersion of Ha (color scale) with the Ha velocity dispersion of 
38 km s -1 overlaid (red contour). Two white lines mark the locations of the 
intersection between x\ and x 2 orbits (contact points). 


timescale on which the two galaxies encounter. Both simula¬ 
tions and observations suggest that the global kinematic dis¬ 
turbances of minor mergers fade within ~ 1 - 2 Gyr after the 
first encounter of two galaxies, and the rotation curves exhibit 
no more severe distortions (e.g. Rubin et al. 1999; Dale et al. 
2001; Kronberger et al. 2006). Thus, NGC 5430 likely captured 
HOLM 569B around a few Gyr ago. However, the merger sys¬ 
tem is not fully settled yet as suggested by the large velocity 
dispersion of the ionized gas, inhomogeneous mass distribution 
at the collision site (the W-R region), and probably the asym¬ 
metric rotation curve at the outer disk as well. 


4 Physical Properties of Molecular Gas and 
Star Formation activity 

4.1 Physical Properties of Molecular Gas 

4.1.1 Line Ratios 

Because of the different excitation conditions associated with 
each molecular line, the line ratio constrains the physical prop¬ 
erties of the molecular gas. In addition to the 12 CO (1-0) to 
13 CO (1-0) ratio (Ri-ol 12 CO (2-1) to 12 CO (1-0) CRiaco) 
and 12 CO (2-1) to 13 CO (2-1) ( R 2 - 1 ) ratios are quoted by us¬ 
ing the measurements reported in the literature. Table 1 shows 
the line intensities and ratios of the galactic center and the W- 
R region of NGC 5430. All values have been measured on a 
scale of 22" (in diameter), centering on the galactic center and 
the W-R region. The scale was chosen by the lowest resolution 
among all measurements, which is the IRAM 30m single dish 
observation in 13 CO (2-1). 

The line ratios of the galactic center of NGC 5430 are similar 
to those of the Galactic star forming regions, starburst nucleus, 


Position [kpc] 

-10 -50 5 10 



Fig. 7. Position-Velocity (PV) diagram of NGC 5430. PV diagram is cut along 
the major axis of NGC 5430 with an position angle of 181.9°. PV diagram 
of 12 CO (1-0) is presented with black contours. PV diagram of Ha is shown 
with color scale and grey contours. The PV diagram of Ha has the same 
spatial and velocity resolution with that of 12 CO (1-0). Blue curve denotes 
the rotation curve traced by Ha created by Epinat et al. (2008). For compar¬ 
ing between the two sides of the galaxies, reserved rotation curve is overlaid 
with pink dots. 


and unresolved starburst galaxies, while the values for the W- 
R region are consistent with the galactic disk clouds. R 1-0 of 
the galactic center and the W-R region of NGC 5430 are ~ 10 
and ~ 6, respectively. The high R 1-0 (10 - 20) is observed 
in the Galactic star forming regions, starburst nucleus, and un¬ 
resolved starburst galaxies, while the low R 1-0 (~ 6) occurs 
in the galactic disk, relatively quiescent molecular clouds (e.g., 
Aalto et al. 1995; Paglione et al. 2001; Simon et al. 2001; Tan et 
al. 2011; Papadopoulos et al. 2012). The high R 1-0 is explained 
by the change in opacity of the molecular gas owing to the high 
temperature, velocity dispersion and column density of active 
star forming region and galaxies. f?i 2 CO also increases with 
increasing star formation activity. The variation stems from the 
fact that the upper energy boundary (^ 16.6 K) and the effective 
(including radiative trapping) critical density (~ 230 cm -3 ) of 
12 CO (2-1) are slightly higher than the typical temperature (~ 
10 K) anddensity (~100 cm -3 ) of molecular clouds. Therefore, 
f?i 2 CO increases from ~ 0.5 in quiescent molecular clouds, to 
> 0.7 in dense and/or warm star forming clouds (e.g., Chiar et 
al. 1994; Oka et al. 1998; Sawada et al. 2001; Koda et al. 2012). 
The line ratio for the W-R region (0.54) is similar to those of 
the quiescent clouds, while a higher ratio (0.7), close to that of 
star forming regions and galaxies, is found for the central re¬ 
gion. Finally, R 2-1 appears to increase from < 10 to 10 - 60 
when moving from quiescent to star forming regions and galax¬ 
ies (e.g., Aalto et al. 1995; Sawada et al. 2001; Papadopoulos et 
al. 2012; Israel et al. 2015). The trend is similar to that of R 1 - 0 . 
R 2-1 is available at the central region of NGC 5430. The mea¬ 
sured value of 12.6 was in the range for star forming regions 
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(a) (b) (c) (d) (e) 

Fig. 5. Velocity fields of NGC 5430. Color scale of velocity fields is shown with a range 2700 to 3100 kms -1 . Contours of velocity fields are plotted in steps 
of 2720, 2770, 2820, 2870, 2920, 2970, 3020, and 3070 km s -1 in all panels. The galactic center and the W-R region are marked with crosses in all panels 
as well, (a) Velocity field of 12 CO (1-0) (color scale and contours), (b) Velocity field of 13 CO (1-0) (color scale and contours), (c) Comparison of velocity fields 
of Ha (color scale and black contours) and 12 CO (white contours), (d) Smoothed velocity field of Ha with velocity and spatial resolutions of 12 CO field, (e) 
Residual velocity of the smoothed velocity field of Hoi (panel d) minus the velocity field of 12 CO field (panel a). 


and galaxies. 

4.1.2 Large Velocity Gradient Calculations 

Gas temperature and volume density can be quantified by com¬ 
paring the observed line ratios with those obtained by using 
one-zone large velocity gradient (LVG) models based on local 
photon trapping and the escape probability method (Goldreich 
& Kwan 1974; Scoville & Solomon 1974). In the LVG cal¬ 
culations, molecular excitation conditions depend on the local 
kinetic temperature (Tk), line opacity (column density of CO 
per unit velocity, Nco/dv ), and volume density (nH 2 ), as low- J 
CO is excited by the collisions with H 2 . The CO-H 2 collisional 
cross-sections reported by Yang et al. (2010) were used in the 
present work. 

The column density of CO per unit velocity 
was assumed by using our observations, and it is 
\og(Nco/dv) [cm _2 (Kkms _1 ) _1 ] — 16.6 — 17.2. Thus, the 
calculation was simplified to consider the variables of Tk and 
nn 2 • Assuming the width of velocity range for giant molecular 
clouds associations to be 40 - 100 km s -1 , and assuming the 
surface density of ~ 1000 M 0 pc -2 and a CO-to-H 2 abundance 
of 8 x 10 -5 , the resultant log(iVco/du) is in the above range. 
We note that the following results of our LVG calculations 
were based on the above assumptions, and therefore should be 
treated with caution. 

Figure 8 shows the results of our LVG calculations. Panel 
(a) and (b) show the results for the central and W-R regions, 
respectively. The results are shown as Tk (y-axis) versus nu 2 
(x-axis). The observed line ratios are shown as contours. The 
dash, solid and dotted-dashed contours denote Ri-o, T 12 CO, 


and R 2 - 1 , respectively. R 2-1 is available for the central region 
only. For each line ratio, two contours are plotted, correspond¬ 
ing to \og(Nco/dv) = 16.6 in thin contour and 17.2 in thick 
contour. Solutions were found in the areas enclosed by all con¬ 
tours; these areas are marked in grey. 


The LVG calculations suggest that the bulk molecular gas is 
not heated by the compact starburst in the W-R region, whereas 
the bulk molecular gas is generally warm in the central re¬ 
gion. Bulk molecular gas associated with the W-R region has 
low temperature of ~ 10 K, similar to that of the galactic disk 
molecular gas. Starburst is ongoing in the W-R region (§4.2); 
thus, the temperature is expected to be high, presumably at or 
above ~ 100 K (e.g., WeiB et al. 2001; Matsushita et al. 2010). 
The low calculated temperature implies that the off-center star- 
burst is spatially small so that the bulk molecular gas seen in 
low-resolution observations (average over the scale of kpc) re¬ 
mains cold. The Tk in the central region is 10 - 40 K, warmer 
than the galactic disk molecular gas (~ 10 K). Although the star 
forming activity in the central region is not as drastic as that in 
the W-R region, the distribution of star forming regions is more 
extended. Thus, the molecular gas is generally warm in the cen¬ 
tral region, yielding a high average Tk over the scale of kpc. 


The LVG calculations show that the density of bulk molec¬ 
ular gas traced by the CO lines is 200 -1600 cm -3 in both re¬ 
gions. The gas density is consistent with the effective critical 
density of the CO lines. 
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Table 1 . Line intensities and ratios of the central (second row) and W-R regions (third row). Measurements of 12 CO and 13 CO (1-0) 
are from this work using CARMA telescope. 12 CO (2-1) and HCN (1-0) are obtained from the IRAM 30-m telescope by Contini et al. 
(1997). 13 CO (2-1) is measured using IRAM 30-m telescope as well by Kruegel et al. (1990). The listed numbers have a resolution of 
22" in diameter. The original resolution of HCN (1-0) is 27". The total flux (in Jansky) within the central 27" has been converted into 
the brightness of a 22" beam based on the assumption that the spatial size of dense gas cannot be larger than that of CO emission. 


12 CO (1-0) 

1 a CO (1-0) 

^ CO (2-1) 

ia CO (2-1) 

HCN (1-0) 

O _ 12 CO (1 — 0) 

„ _ 12 CO (2—1) 

— 12 CO (2—1) 

„ _ HCN (1-0) 

dense 12 CO (1-0) 

[Kkms -1 ] 

[Kkms -1 ] 

[Kkms -1 ] 

[Kkms -1 ] 

[Kkms -1 ] 

1 0 13 CO(1-0) 

12CO 12 CO (l 0) 

2 1 13 CO(2-l) 

117.4 ± 2.1 

12.7 ± 0.1 

81.6 ± 0.5 

6.5 

2.1 ± 0.3 

10.3 ± 0.1 

0.70 ± 0.01 

12.6 

0.018 ± 0.001 

11.4 ± 0.3 

1.9 ± 0.1 

6.1 ± 0.5 


< 0.8 

6.0 ± 0.1 

0.54 ± 0.01 


< 0.070 



volume density n(H z ) [cm' 3 ] 


(a) 



volume density n(H z ) [cm' 3 ] 

(b) 

Fig. 8. Result of LVG calculations in (a) the central region and (b) the W- 
R region. The dashed, solid and dotted-dashed contours represent the 
ratios of 12 CO(1-0)/ 13 CO(1-0) (i?i_ 0 ), 12 CO(2-1)/ 12 CO(1-0) (i?i 2 co), 
and 12 CO(2-1)/ 13 CO(2-1) ( R 2 -i ), respectively. Ratio of R 2 -1 is available 
for the central region only. For each ratio, two contours are drew, repre¬ 
senting the \og(Nco/dv ) = 16.6 (thin contour) and 17.2 (thick contour) 
cm -2 (Kkms -1 ) -1 , respectively. Solutions of LVG calculations are high¬ 
lighted with grey shadows. 


4.2 Star Formation Activity 

4.2.1 Star Formation Rate 

History of star formation can be quantified by using star forma¬ 
tion rates (SFRs) tracing different timescales. In this work, we 
used the SFRs estimated from the infrared (SFR(IR)), radio con¬ 
tinuum (Figure 9(a), SFR(RC)), and Ka (Figure 9(b) and 9(c), 
SFR(Ha)) to constrain the history of star formation. The three 
wavelengths respectively trace the average SFRs during the last 
~ 100 Myr, several 10s of Myr, and several Myr, respectively 


(Kennicutt 1998a; Murphy et al. 2011; Calzetti 2013). All SFRs 
in this work were estimated based on the Kroupa initial mass 
function (IMF). 

Infrared observations suggest that the average SFR during 
the last ~ 100 Myr is significantly higher than that of the Milky 
Way. Infrared emission originates from the dust heated by stars 
younger than 10 8 yr; therefore, SFR(IR) traces a relatively long 
history of star formation. The luminosity-to-SFR relation, cali¬ 
brated by Calzetti (2013), is: 


SFR 

[M e yr-i] 


= 2.84 x 10“ 44 


Lm 

[ergs- 1 ] ’ 


(2) 


where Lm is the infrared bolometric luminosity in the 5 to 1000 
pm range. The Lm of NGC 5430 is 10 10 ' 88 L 0 as obtained 
from IRAS, resulting in the SFR(IR) of ^ 8 M 0 yr -1 , ~ 8 
times higher than that of our galaxy. 

The value of the SFR(RC) suggests that an instantaneous 
starburst had occurred ~ 10 Myr ago. Radio continuum at 1.4 
GHz is mainly the non-thermal synchrotron radiation (~ 90%) 
produced by relativistic electrons ( S ). These electrons are ac¬ 
celerated by the supernovae of massive stars. Therefore, ra¬ 
dio continuum traces the star formation history on the timescale 
corresponding to the lifetime of the lowest massive stars (8 M 0 ) 
that could become supernovae, which is ~ 40 Myr. Condon et 
al. (2002) and Heesen et al. (2014) suggested a derivation of 
SFR based on the radio continuum: 


SFR 

[Moyr- 1 ] 


= 0.75 x 10“ 21 


£l.4GHz 

[WHz -1 ] ’ 


(3) 


where Li. 4 gh z is the luminosity in at 1.4 GHz. A small fraction 
of radio continuum originates from free-free emission (Sth). 
We use the following relation: 


to subtract the thermal emission from the total flux, where a 
is the non-thermal spectral index and v is the observed fre¬ 
quency. We adopt a typical value of a « 0.8 (Condon et al. 
2002). The total flux at 1.4 GHz measured by the VLA was 65.9 
mJy (Condon et al. 1998). After subtracting S t h from Si.4gh z , 
f i. 4 GHz becomes 1.3 x 10 22 W Hz -1 . The corresponding SFR 
is ~ 10 M 0 yr -1 , one order of magnitude higher than that of 
the Milky Way. Among the global SFR(RC), ~ 1 M 0 yr -1 is 
contributed by the W-R region. The instantaneous starburst may 
be a result of the satellite collision. Simulations predict that the 
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highest SFR occurs during the final merging phase of two galax¬ 
ies and during the post-merging phase (Cox et al. 2008; Lotz et 
al. 2010). 

The derived SFR(RC) is a lower limit owing to the mass loss 
of massive stars in the W-R region. Keel (1982) argued that the 
mass loss of massive stars in the W-R region is too extreme, so 
that only ~ 4% of the massive stars could become supernovae 
and contribute to the radio continuum. If true, the SFR(RC) is 
considerably underestimated. If we naively scale the observed 
SFR(RC) by the fraction of stars that could not become super¬ 
novae, the real global SFR during the past ~ 10 Myr will be 
around 10 (the observed global SFR(RC)) + 1 (the observed 
SFR(RC) in the W-R region)/0.04 = 35 M 0 yr“\ and ~ 25 
M 0 yr -1 will be contributed by the W-R region. 

In contrast to the past, the recent SFR (few Myr) indicated 
by the total luminosity of the ionized gas is not particularly 
dramatic, comparable to that of the Milky Way. The Ha from 
the ionized gas traces the most recent star formation owing to 
the short lifetime of massive stars. Figure 9(b) shows the Ha 
image obtained from the Observatoire Haute-Provence 1.92-m- 
diameter telescope as part of the kinematical 3D Gassendi Ha 
survey of spirals survey (GHASP). The Ha luminosity (Lh«) 
suggested by the survey data is 2.4 x 10 41 erg s -1 (Epinat et 
al. 2008). SFR(Ha) is derived from the luminosity-SFR cali¬ 
bration of Calzetti (2013): 


SFR 

[M 0 yr-i] 


= 5.5 x 10“ 42 


Lhu 

[ergs- 1 ]' 


(5) 


The derived global SFR (Ha) is 1.3 M 0 yr -1 . This value sug- 
gests that the recent global SFR is comparable to that of the 
Milky Way. 

The Galactic level SFR is confirmed by adding up the lo¬ 
cal extinction-corrected SFR of individual HII regions iden¬ 
tified by other independent observations. Figure 9(c) shows 
the HII regions identified by Briere et al. (2012) by using the 
imaging Fourier transform spectrograph SpIOMM on the 1.6- 
m-diameter Ritchey-Chretien telescopes. The total SFR in the 
HII regions is 2.4 M© yr -1 , including ~ 0.6 M© yr -1 from 
the central region and 1 M© yr -1 from the W-R region. The 
values are higher than the SFRs suggested by GHASP obser¬ 
vations, probably owing to the uncertainty of the indirect flux 
calibration in the GHASP survey (see Epinat et al. 2008 for 
the details). In spite of the difference, both global SFR(Ha) 
are significantly lower than the global SFR(IR) and SFR(RC). 
Moreover, the SFR in the W-R region decreases by more than 
10 times between the actual SFR(RC) and SFR(Ha). 

NGC 5430 is experiencing the end of the recent W-R star- 
burst. Evolutionary stellar synthesis studies suggest that the age 
of the W-R region in NGC 5430 is > 5.3 Myr (Fernandes et al. 
2004). It is at the relatively late stage because W-R starbursts 
typically fade out within 10 Myr (Schaerer & Vacca 1998). Our 
derived SFRs also suggest that the star formation activity had 


peaked around ~ 35 M© yr -1 , and is rapidly decreasing to 
the Milky Way value (1-2 M© yr -1 ). Moreover, theoreti¬ 
cal models suggest that the ionized flux generated by massive 
stars decreases by 1 - 2 orders at the end of a W-R starburst 
(Schaerer & Vacca 1998). The difference between the SFRs of 
actual SFR(RC) and SFR(Ha) from the W-R region is consis¬ 
tent with this theoretical prediction. 

Finally, if the W-R episode occurred immediately after the 
satellite collision, the collision time should be between ~ 5.3 
and 10 Myr ago because the W-R period (the time during which 
the ionized flux decreases to the pre-starburst level or the dura¬ 
tion of the appearance of the W-R stars) theoretically does not 
exceed 10 Myr (Schaerer & Vacca 1998). 

4.2.2 Relation Between the Current Star Formation 
Activity and Molecular Gas 

In this section, we compare star formation activity and molec¬ 
ular gas. Since the Ha flux is not spatially calibrated on the 
image of the GHASP survey (Figure 9(b)) due to the strategy 
of the survey, we used the calibrated flux of HII regions from 
Briere et al. (2012) and the 12 CO flux of the beam at the same 
positions for the analyses in this section. Typical size of HII 
regions is 1 - 2.5 kpc 2 , comparable to our beam area of ~ 1.5 
kpc 2 , but we suggest that the results still have to be treated with 
care. 

We found that the 12 CO distribution of NGC 5430 is corre¬ 
lated with star formation rate surface density. Figure 9(d) com¬ 
pares the HII regions of Briere et al. (2012) (circles) and our 
12 CO image (color scale). Physical areas of HII regions are in¬ 
dicated by the circles sizes. To fairly compare the star formation 
activity, we normalized the SFRs of HII regions by the their ar¬ 
eas, to obtain the quantity that is known as the star formation 
rate surface density (Esfr), measured in the units of M© yr -1 
kpc -2 . The Fsfr of the HII regions is indicated with colored 
circles. Cyan, yellow and magenta colors denote the Esfr be¬ 
low 0.005 M© yr -1 kpc -2 , in the 0.005 - 0.05 M© yr -1 kpc -2 
range, and above 0.05 M© yr -1 kpc -2 , respectively. The 12 CO 
emissions are mostly associated with moderate and high Esfr 
(yellow and magenta circles). The missing 12 CO in the low 
Esfr regions can be explained by insufficient sensitivity, by 
assuming the ratio of Esfr to 12 CO flux to be the same in the 
moderate- and Iow-Esfr HII regions (in other words, the same 
star formation efficiency). 

Current global star formation efficiency (SFE) of NGC 5430 
is consistent with those of other disk galaxies. The current SFE 
is quantified as SFR(Ha)/MH 2 , where Mh 2 is the molecular gas 
traced by the 12 CO. The total SFR(Ha) in the 12 CO region is ~ 

2.2 M© yr -1 obtained by using the measurements of Briere et 
al. (2012). The total Mh 2 in our map is ~ 4.7 x 10 9 M©, for 
the total flux of 5.7 x 10 4 Jy beam -1 km s -1 and the CO-to-H 2 
conversion factor of 2 x 10 2 ° cm -2 (Kkm s -1 ) -1 . The global 





12 


Publications of the Astronomical Society of Japan, (2014), VoI. 00, No. 0 



1010 0.00029 0.00088 0.00187 0.00324 0.0C ) 3 14 32 56 88 126 171 224 284 3‘ ) 3 14 32 56 88 126 171 224 284 3‘ 1 0.8 3.2 7.2 13 20 29 39 51 65 8 


(a) (b) (c) (d) 

Fig. 9. Star formation tracers of NGC 5430. Multi-wavelength intensity maps are shown in color scale and contours, where contours are plotted with steps of 
10, 15, 40, and 70 % of the maximum intensity. Resolution of the observations is shown in the lower-right corner of each panel, (a) Radio continuum in 20 cm 
(1.4 GHz) observed by VLA (Condon et al. 1998). (b) Ha image from the GHASP survey (Epinat et al. 2008). (c) Hll regions (circles) identified by Briere et 
al. (2012) overlaid on panel (b). Size of circles denote the physical size of Hll regions, (d) Hll regions (circles) of Briere et al. (2012) overlaid on 12 CO image. 
The size and color of circles indicate the physical size and star formation rate density (E S fr) of Hll regions, respectively. The colors denote, magenta: £ S fr 
> 0.05 M Q yr -1 kpc -2 , yellow: 0.005 < £ S fr < 0.05 M 0 yr -1 kpc -2 , and cyan: £ S fr < 0.005 M 0 yr -1 kpc -2 . 


SFE of the galaxy is therefore ~ 4.6 x 10“ 10 yr -1 . The value 
is consistent with those for the nearby disk galaxies (Kennicutt 
1998b; Bigiel et al. 2008; Huang & Kauffmann 2015). 

The local SFE exhibits large regional variation. The SFE at 
the galactic center is ~ 9 x 10 10 yr 1 . The SFE is calculated 
by using the total SFR(Ha) of the central two HII regions from 
Briere et al. (2012) and total 12 CO flux of the beams associated 
with them. We found that the SFE is slightly lower compared 
with other circumunuclear regions, for which the SFEs are typ¬ 
ically on the order of 10 -9 yr -1 (e.g., Kennicutt 1998b). 

The SFE of the W-R region is 10 times higher than the global 
SFE and the SFE at the galactic center, and is ~ 4.5 x 10“ 9 
yr -1 , as obtained by using the SFR(Ha) of the giant HII region 
at the W-R region and the 12 CO flux of the beam associated with 
it. The large beam of CO observations may include molecular 
gas irrelevant to the starburst. Therefore, the derived SFE is a 
lower limit. 

4.2.3 The Externally Triggered Starburst at the Bar End 

Watanabe et al. (2011) suggested that the enhancement of star 
formation efficiency at the bar ends is owing to the effective 
cloud-cloud collisions triggering the star formation. This ques¬ 
tions whether HOLM 569B is indispensably required to explain 
the W-R starburst at the eastern bar end. W-R stars do appear at 
the bar end of the Milky Way and are more concentrated towards 
the spiral arms, suggesting a possible promotion of massive star 
formation via galactic structures (e.g., Davies et al. 2012). The 
western bar end is lacking unusual star formation. Because the 
W-R starburst is transient, it is possible that it had taken place 


and terminated at the opposite side end. In other words, the 
two bar sides do not evolve simultaneously (Keel 1982). This is 
certainly a possible alternative. 

Nonetheless, the W-R region exhibits several star formation 
features that do not favor an internal origin. The surface den¬ 
sity of the W-R stars in the W-R region of NGC 5430 is more 
than 100 times higher than that of the Milky Way (Massey 
2003; Mauerhan et al. 2011; Davies et al. 2012; de la Fuente 
et al. 2013). Moreover, the HII region associated with the W- 
R region deviates from the global HII regions in terms of the 
luminosity function, which is similar to our Milky Way (Keel 
1982). Such a distortion of the HII regions luminosity func¬ 
tion is not observed in the internally triggered extreme starburst 
environments such as circumnuclear starbursts (e.g., Kennicutt 
et al. 1989; Feinstein 1997; Alonso-Herrero & Knapen 2001). 
Finally, the star forming region in the W-R region is particularly 
young (~ 5 Myr) compared with the rest of the bar area, which 
is 10 - 12 Myr old. All these features suggest that the W-R re¬ 
gion has been likely formed by the collision with HOLM 569B. 

4.2.4 Dense Gas Content and Future Star Formation 
Activity 

Dense gas content (> 10 4 cm -3 ) of molecular gas can be es¬ 
timated in terms of dense gas fraction.Table 1 lists the inten¬ 
sity of dense gas HCN (1-0) and the ratio of HCN/ 12 CO (1-0) 
(Odense) in the central and W-R regions, Rdense increases from 
0.01 - 0.06 in non-starburst nuclei to 0.06 - 0.2 in starburst nu¬ 
clei (Matsushita et al. 2010). i^dense in the central region of 
NGC 5430 is within the range of non-starburst nuclei, 0.018. 
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The value for the W-R region remains undetermined, < 0.070. 
Odense is often used to infer the fraction of dense gas (/dense) 
because the critical density of HCN (1-0) is as high as 10 4-5 
cm -3 . Together with the adopted Xco and the HCN-to-H 2 
conversion factor (X H cn) of 1.3 x 10 21 cm -2 (K km s -1 ) -1 
(Solomon et al. 1992), Flense = 0.018 corresponds to /dense as 
low as ~ 12 %. In the W-R region, /dense is nearly four times 
higher than that in the central region (46%), by using the upper 
limit of 0.070 on R d ense* 

The dense gas fraction in the central region of NGC 5430 
is lower than that in the circumnuclear regions of other strong- 
barred galaxies (~ 50 - 100%, e.g., Kohno et al. 1999; Pan 
et al. 2013), implying that even though the strong bar is effi¬ 
ciently transporting molecular gas to the center, the majority of 
the gas does not become dense gas. The low density gas frac¬ 
tion explains the above-mentioned low SFE of the central region 
compared with other circumnuclear regions. According to the 
relation of HCN luminosity to SFR that was suggested by Gao 
& Solomon (2004), the future SFR formed by this dense gas 
will be ~ 0.4 M 0 yr -1 , comparable to the current SFR in the 
central region indicated by Ha. 

The future SFR, arising from the dense gas in the W-R re¬ 
gion, will be at most 0.15 M 0 yr -1 , nearly 10 times lower than 
the current SFR in the W-R region, which is ~ 1 M 0 yr -1 . 
This is consistent with the previous discussion suggesting that 
the star formation activity in the W-R region is decaying and 
will continue to decrease in the future. 

5 Initial Mass Ratio 

We roughly estimate the initial mass ratio of NGC 5430 and 
its companion based on the observed morphology, velocity, and 
star formation properties. Major mergers can destroy galaxy 
disks and produce remnants with properties similar to those of 
observed giant ellipticals (but some studies have reported the 
possibility of formation of a disk galaxy in a major merger (e.g., 
Springel & Hernquist 2005). We rule out this possibility since 
the candidate of companion has been suggested to be embedded 
in NGC 5430 (Keel 1982). Intermediate-mass mergers (4:1 - 
7:1) tend to produce both ellipticals and hybrid systems, which 
are disk-like galaxies supported by velocity dispersion rather 
than rotation (Bournaud et al. 2004; Bournaud et al. 2005). The 
value of rotation velocity to velocity dispersion of these hy¬ 
brid systems (V ro t /cr) is small, 0 < V YO t/<J < 2 (Jog & Chitre 
2002; Naab & Burkert 2003; Bournaud et al. 2004; Bournaud et 
al. 2005). The maximum Ha velocity dispersion of the galaxy 
is around 70 km s -1 , and the deprojected Ha rotational velocity 
in the galactic disk is ~ 300 km s -1 (corrected with an inclina¬ 
tion of 32° from Epinat et al. 2008), resulting to a Vrot/cr of 
~ 4. Hence, NGC 5430 is strongly rotationally-supported as 
a regular disk galaxy. Therefore the initial mass ratio between 


NGC 5430 and its satellite is likely larger than 7:1. Simulations 
show that, for a primary galaxy with a stellar mass of NGC 5430 
value 5 , enhanced star formation only occurs in pairs with mass 
ratio < 20:1 (Cox et al. 2008; Lotz et al. 2010). Therefore the 
initial mass ratio of NGC 5430 is likely between 7:1-20:1, and 
probably around 10:1 to generate a star formation timescale (in¬ 
verse of SFE) of < 2 Gyr at the initial coalescence phase (Cox 
et al. 2008). 

6 Summary 

In this paper, we constrain the minor-merging and star forma¬ 
tion history of the three-spiral Wolf-Rayet (W-R) galaxy NGC 
5430 and its satellite HOLM 569B embedding in the eastern bar 
of NGC 5430. 

Molecular gas tracers 12 CO (1-0) and 13 CO (1-0) were ob¬ 
served by CARMA with a resolution of ~ 1.3 kpc (§2). The 
observed area is ~ 16 kpc. The main results are (§3): 

• 12 CO (1-0) is confidently detected in the central region, en¬ 
tire bar, and an asymmetric spiral arm (the third spiral) of 
NGC 5430. Solid detections of 13 CO (1-0) show at the cen¬ 
tral region and the eastern bar with W-R region. 

• The molecular bar shows as offset ridges as seen in typical 
barred galaxies. The offset ridges of the two sides of the bar 
are symmetric around the galactic central region, but the east¬ 
ern bar are stronger in both CO lines. 

• NGC 5430 has a strong bar suggested by the high central 
concentration of molecular gas and the high ellipticity of the 
bar. 

• Velocity field of 12 CO (1-0) matches the velocity field and 
the Position-Velocity diagrams of Ha, which was observed 
with higher spatial and velocity resolutions. 

• Rotation curve of Ha shows no obvious distortion by galax¬ 
ies merger. However, inhomogeneous mass distribution is 
suggested by a small bump of rotation curve around the W-R 
region, implying that the merger system has nearly, but not 
fully, stablized yet. The outer rotation curve are more asym¬ 
metric between the two sides of the galaxy due to either the 
past interaction with HOLM 569B or/and ongoing interaction 
with members in the galaxy group. 

• Since the minor merger systems kinematicaly stablize after 1 
- 2 Gyr after the first encounter, NGC 5430 likely captures 
HOLM 569B on this timescale. 

Physical properties of molecular gas at the central region 
and the W-R region are inferred by the Large Velocity Gradient 
calculations (LVG) and the line ratios of 12 CO (1-0)/ 13 CO (1- 
0) CRi_o), 12 CO (2-l)/ 12 CO (1-0) (ifoco), and 12 CO (2- 
l)/ 13 CO (2-1) (R 2 - 1 ) (§4). The later two ratios are cited from 

5 Stellar mass of NGC 5430 is ~ 5.3 x 10 10 M o from the Spitzer S 4 G project 
(Sheth et al. 2010). 
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literature. The ratio of R 2-1 is only available in the central re¬ 
gion. The main results are: 

• The line ratios of the central region are R 1-0 « 10, f?i 2 CO ~ 
0.7, and R 2-1 ~ 12.6. The ratios are in the range of starburst 
regions, but close to boundaries between starburst and the 
relatively quiescent molecular regions (e.g., the Galactic disk 
molecular gas). The ratios of the W-R region are consistent 
to the galactic disk molecular gas, with R 1-0 and i?i 2 CO and 
~ 6 and 0.54, respectively. 

• Temperature of the central region derived from the LVG cal¬ 
culations is 10 - 40 K, warmer than the W-R region, which 
is 10 K. These results imply that the W-R region is small 
and compact, thus the bulk molecular gas measured by low- 
resolution observations includes significant amount of low- 
temperature molecular gas irrelevant to the W-R starburst. In 
contrast, in spite of a lower degree of star formation, star for¬ 
mation is more spatially extended in the central region, lead¬ 
ing to a higher derived temperature. 

• Gas density is 200 - 1600 cm -3 in both regions, consistent 
to the critical density of CO lines. 

Star formation history is constrained using star formation 
rates (SFRs) that trace different timescales, including infrared 
(SFR(IR)) for the average of past 100 Myr, radio continuum 
(SFR(RC)) for the past a few times 10 Myr, and Ha (SFR(Ha)) 
for the past a few times 1 Myr (§4). The main results are: 

• The observed global SFR(IR), SFR(RC) and SFR(Ha) are ~ 
8, 10, and 2.4 M© yr -1 , respectively. The observed SFR(RC) 
might be an lower limit due to the mass loss of massive stars. 
The real SFR(RC) is probably as high as ~ 35 M© yr -1 . 
These SFRs suggest that there is an instantaneous starburst 
around a few times 10 Myr ago. Together with the previous 
stellar synthesis study of the W-R region and the theoretical 
model of W-R starburst, we suggest that the starburst started 
around 5.3-10 Myr ago, likely triggered by the collision of 
the satellite. 

• We found that the current star formation efficiency of the cen¬ 
tral region of NGC 5430 is low compared to other galaxies. 
Moreover, the current dense gas fraction of the central region 
inferred by the ratio of HCN/ 12 CO (1-0) is also low, imply¬ 
ing a low SFR in the future as well. 

• Current SFE of the W-R region is 10 times higher than the 
central region. However, SFR of the W-R region is rapidly 
decreasing by more than 10 times from the collision time to 
the present time. The future SFR implied by the amout of 
dense gas suggests that the SFR of the W-R region will con¬ 
tinue to decrease. 

We discussed the possible initial mass ratio of NGC 5430 
and its satellite to be 7:1 - 20:1 based on the observed morphol¬ 
ogy, velocity, and star formation properties (§5). 

Finally we note that theoretical simulation is still necessary 


to constrain the merger history of NGC 5430 system. Each 
systems may have unique satellite orbits that would signifi¬ 
cantly affect their current morphology and star formation ac¬ 
tivity. Models of NGC 5430 and the direct comparison of the 
multi-wavelength observations will yield important constraints 
to the this system. Future high resolution HI observations would 
be valuable to reveal large scale dynamic structures beyond the 
molecular and stellar disk and provide a knowledge of the pos¬ 
sible interaction orbits. Wide-field CO observations are also in¬ 
dispensable to reveal the gas distribution and gas supply to star 
formation in detail, as well as the possible influence of galaxies 
collisions towards the central region, e.g, the asymmetric SFE 
at the contact point of the circumnuclear ring. 
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